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This publication reports research involving pesticides. It does not contain recommendations for 
their use, nor does it imply that the uses discussed here have been registered. All uses of 
pesticides must be registered by appropriate State and/or Federal agencies before they can 
be recommended. CAUTION: Pesticides can be injurious to humans, domestic animals, 
desirable plants, and fish or other wildlife—if they are not handled or applied properly. Use all 
pesticides selectively and carefully. Follow recommended practices for the disposal of surplus 
pesticides and pesticide containers. The Forest Service, United States Department of Agriculture, 
has developed this information for the guidance of its employees, its contractors, and its 
cooperating Federal and State agencies, and is not responsible for the interpretation or use of 
this information by anyone except its own employees. The use of trade, firm, or corporation 
names in this publication is for the information and convenience of the reader, and does not 
constitute an endorsement by the Department of any product or service to the exclusion of 
others that may be suitable. 


The United States Department of Agriculture (USDA) prohibits discrimination in its programs 
on the basis of race, color, national origin, sex, religion, age, disability, political beliefs, and 
marital or familial status. (Not all prohibited bases apply to all programs.) Persons with 
disabilities who require alternative means of communication of program information (braille, 
large print, audiotape, etc.) should contact the USDA Office of Communications at (202) 720- 
2791. To file a complaint, write the Secretary of Agriculture, U.S. Department of Agriculture, 
Washington, D.C. 20250, or call 1 (800) 245-6340 (voice) or (202) 720-1127 (TDD). USDA is 
an equal employment opportunity employer. 
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Introduction 





A typical pesticide spray pattern over a forest in uniform 
terrain is depicted in Figure 1, where A is the target area 
and the lines represent the flight path of the pesticide 
application aircraft. The aircraft creates trailing vortices, 
which greatly influence the initial state of the droplets. We 
assume here that these vortices have decayed. As the 
pesticide aerosol evaporates, the changing size of the 
droplets and their changing aerodynamics can be 
characterized. Depending on the dispensing nozzle and 
initial vortex mixing, there will be a characteristic droplet 
size distribution over a given initial volume; the larger 
droplets will separate from the smaller ones, 
spatiotemporally altering the droplet size distribution. 
Consider a spatiotemporally varying droplet distribution 

c (r, x, Z, t), where rdenotes the radius of the droplets, x is 
the streamwise (with the wind) space variable, zis the 
vertical space variable, and tis time. The problem may be 
generalized to include a third space variable y and 


corresponding velocity V, to admit the situation where the 
spray pattern is not perpendicular to the wind direction. 
Figure 2 shows the definition of the variables relating to 
height A of the initial droplet distribution, tree height H, and 
5. There are inherently two streamwise length scales in this 
problem. One scale, 5, is associated with the distance 
between each pass of the aircraft through the target region. 
The other scale, the streamwise length 6, of the entire 
target region, is associated more with the off-target drift of 
the smallest droplets and the pesticide vapor. For an 
appropriate choice of spray nozzle and wind conditions, the 
greater proportion of the droplets are immobilized in the 
vegetation in a very restricted region around the area A in 
Figure 3, but a potentially significant amount of off-target 
pesticide can be windborne for great distances. The 
simulation of the larger scale problem is of importance in 
risk assessment and spray optimization; in this paper only 
the smaller scale problem will be addressed directly. 





Figure 1—Aerial pesticide application pattern over forest area “A,” showing definition of distances 4, é. 
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Figure 2—Initial distribution ® of pesticide showing periodic pattern and def 





Figure 3—Off-target pesticide drift at three different times. 
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Physical Parameters 
and Assumptions 





Saturation and Temperature—The results of a pesticide 
application can be affected by the ambient saturation and 
temperature, which in these models are taken to be 
constant through the vertical range, assuming that the 
spray is applied near the treetops. In pesticide evaporation 
models, a primary constant controlling the dynamics is 


ae ee SH et 
ie be flee 
KRT? De.(T) 


(of units [r?/f]) where S is saturation, L is latent heat of 
vaporization, p, is liquid density, K is air thermal 
conductivity, Ry is the individual gas constant for the vapor, 
T is temperature, Dis the diffusivity of the vapor, and e, is 
the enthalpy. These parameters may be estimated from 
ambient conditions and thermodynamic relations for the 
pesticide evaporating in air. If the saturation S>1, we may 
expect condensation to be important; if S<1, we may expect 
evaporation to be important. 


Dilution—The droplets are assumed dilute enough in the 
air that air velocities are unaffected by the droplet 
velocities, and droplets rarely collide. The pesticide vapor is 
also assumed dilute, such that local vapor concentrations 
do not affect the value of a. 


Turbulence and Wind Velocities—Turbulent dispersion is 
neglected through dropping the second-order derivatives in 
the dynamics equations (that is, taking the Sherwood 
numbers to be infinite); the dispersion considered arises 
purely from the average turbulent shearing profile of the 
streamwise wind, and from droplet dynamics including 
sedimentation and evaporation. This is least accurate for 
higher average wind speeds and the smallest droplets or 
vapor. For the purpose of illustration, the vertical 
dependence of the average wind velocity is assumed to be 
reasonably fitted by an elementary function, here assumed 
to be exponentially changing to a constant as height 
increases: 


°° 


U(z) = U_(1- e 7) (1) 


Droplet Deposition onto Foliage—Although the 
deposition of droplets onto foliage is radius-dependent over 
the range of interest, it is assumed here to be dependent 
only on the vertical spatial coordinate. The rate of loss of 
droplet concentration to deposition is assumed to have a 
linear dependence on droplet concentration. The deposition 
coefficient is assumed to vary directly with foliage density. 
We assume a simple form 


k(z)='k,e? (2) 


where k, has units [7/4]. 


Droplet Evaporation—The derivative d//dt describes the 
change in droplet radius with time due to evaporation (or 
condensation), and is approximately given by [1, 2] 
(ignoring the ventilation factor), 


For large radius droplets this can be simplified, disregarding 
solution and curvature effects, to yield the usual 
approximation 


To OL. (3) 


Droplet Velocity—The vertical velocity will be considered 
linear in the droplet radius dependence 


Ol = ob fe 


where Q, has units [z/(rt)]- More accurately, for water 
droplets between 10 and 40 um, there is a quadratic 
dependence; the dependence is linear from about 40 to 600 
um, and then proportional to r”? for larger droplets. These 
nonlinear deviations are assumed negligible. 


Transport Equations 





The transport model for the droplet concentration c (r, x, Z, t), of units (droplets/unit volume) is: 
oc oq , +o(F -c)-50,2 ap # (5a) 


SA EL iC PI 
eS ee Pie Beene 


Here D; are the droplet turbulent dispersion parameters, and q is the concentration of deposited droplets within the vegetative 
canopy. 


The vapor concentration, C(x, z, t) (of units [mass/unit volume]), from evaporated pesticide satisfies 


Pet v-Vo, 1,20) le XZt)+1(% 2.0) Tay Cea ae 
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where |,, |,, |, are source contributions from free droplets, deposited droplets, and infinitesimal droplets (from the edge of the 
droplet spectrum) respectively. These source terms are here assumed to satisfy 


= 4na0| ¢r’ Clr’, Nez) ie (6a) 


ot 
al, (6b) 
— = 4na0| § ror, x,2t\0ms 

a (6c) 


r>0 


xe =i) Cle x 2; | 


The parameter o is the conversion factor of mass of vapor created per unit volume of liquid evaporated, and is specific to the 
pesticide and temperature. Equation (6b) reflects the simplifying assumption that deposited droplets have surface area as if 
they are spherical and entirely exposed to the air. 


Under the given assumptions the droplet equation (5a) is 


dc eV OC oc , a1 O - 
seetat (1-e* 5 + Qor = + koe 745 (%e)=0 (7) 


The total droplet number loss from the target area defined by c(r, x, Z) €[0, -]x[0, 5]x[0, o-] is characterized by the change of 


ayaa ie feel c(r,X,Z)dz ax dr; (8) 


the loss fluxes can be calculated from integration of (7) over c(r, x, Z) (not including diffusive fluxes): 


“ = J Fl 5 (| 02 dr -| at Ue 02 dr + | a 3(Q2C| 9 ax dr, (9) 





that is, from the edge of the droplet spectrum, the horizontal drift, and sedimentation to the soil, respectively. 


All the vapor can be considered to be mass lost “off target,” calculated as a total from (6a, 6b, 6c) 


Vanes (ial fees srldxdzect (10a) 


because it is assumed not to be deposited in the foliage. Similarly to the terms in (9), other mechanisms whereby pesticide 
mass can be off target are from droplets drifting out of the target region, whether horizontal loss 


4n fg oe 
Dor = a Pp| 0 | 0 i} r°c(r,6, Z,t)dr dz dt (10b) 
or vertical loss 
or aul } al al oer, x, 0, Hamaxeat (10c) 
3 p 0 0 0 * 


where P,, is the liquid pesticide density, and diffusive flux has been excluded. 


In order to set the problem in perspective, the second-order equation (5a) is cast also in its dimensionless form: 


OC OC. a OC me cand [C OMOCl oO. ~OC 
ee tes ah eG Cog taco <4 oko ge nee gid (11) 
oo ala 5+ aT on 72 (fre Ey: yeaa ie Taoh 


where d,=D,/8k, is the inverse effective horizontal Sherwood number, d,=D,/H’k,, is the inverse effective vertical Sherwood 
number, t=t/k,, €=x/8, C=z/H, p=rr_,, where r,_,, is the approximate maximum size of droplet (always limited by breakup 
instability), and the Stanton numbers are S,=5k/U_ (horizontal), S.=Hk,/(Q,/,,,,.) (vertical), and S, (evaporative). Other 
dimensionless numbers of importance are 8/H and e=Q,/U.. 


There are various parameter limits of physical interest, namely Sherwood numbers Sh-30, ~ and velocity ratios e0, ». The 
hyperbolic problems considered here all satisfy Sh--, appropriate for ballistic and almost ballistic droplets. Also, we 
consider A/H=2, illustrating a case where the aircraft has flown fairly close to the canopy height. 


Boundary and initial conditions for the hyperbolic equation (7) are c-0 for x=teo, and c (r, x, z, 0) = o(r, x, Z). Specifically, the 
initial condition can be modeled as (z-A)? ny (x-j5)? 


_—br- 
g(r,x,z) = Ar’e 


40, y @ 1 
j=0 
but in the examples below we focus on the cloud of pesticide from only one pass of the aircraft, given in the form 
2 2 
br) — 
g(r,x,z) = Arve : Be 





(12) 


or, undimensionalized, 
(6-2/6)? _(¢-A/H)? 
ee 4o- 4o- (13) 


0(p,6,6) = Ap’e 





This is in accordance with reducing the number of possible solutions to the optimization problem, disregarding diffusive 
transport, and is consistent with keeping the choice of U_ to a relatively small value for pilot safety and reduction of off-target 
drift and off-target deposition. We take 6,, ©,, k, to be fixed; but more generally the adjustable parameters for optimization 
include n, b, A, A, 5, U_, a, and k,. 


“Cost” Functions 





In order to specify an optimization problem, one must characterize a notion of “how close to optimum” a particular solution is. 
With this in mind, introduce a function P=P(z), which is the average (measured or theoretical) distribution of pests in the 
foliage. Some pests are concentrated in low bushes, while others are concentrated in the upper foliage of trees. 


The function P can then be compared to an average dosage measure D(z) by a “mortality” function M(z)=aP(z), simply 
proportional to P and indicating the desired pesticide mass dosage according to pest distribution and pesticide efficacy. Dis 
calculated relative to a particular solution to the transport equations: 


4n a a 
D(z)e a Da % | | =rq(rexyz, Hardxdt: 


The comparison 
C.(D) = J ,(M(z) -D(z)) dz (14) 


then can be minimized with respect to whatever free parameters appear in D(z). Typically, one must solve the optimization 
problem numerically. 
An important consideration is that of the off-target pesticide. Two major off-target effects are possible: 

e@ If the pest is only in the vegetation, then pesticide reaching the soil is off target. 

@ Pesticide drifting horizontally may reach locations outside the prescribed application area. 
Whether or not the “on-target” cost function C, is minimized, the off-target drift may be large or small, and might itself need to 
be minimized, particularly if the pesticide is noxious. If the off-target drift is to be minimized, then for a particular solution, we 


calculate off-target dosage D,, as in (10), which might include both adsorbed and unadsorbed pesticide, perhaps including soil 
uptake and vapor production, and minimize 


Cy= Dept Dor + Vor (15) 


To simultaneously minimize on-target dosage error and off-target drift, one may define a “penalty function” L and “penalty” y 
such that 


On yee 


The penalty y can be adjusted to change the relative weight of importance of the two cost functions C, and C.,. 


Examples 





In this section we consider some aspects of the pesticide 
application problem. (Analytical results are listed in the 
appendix.) Figures 4 and 5 show the evolution of an initially 
symmetric droplet profile, including free and deposited 
droplets. In Figure 4, the horizontal Stanton number S, is 
large, while in Figure 5 it is small, showing more of an effect 
of the shearing wind. If we look at the total pesticide mass 
deposited as only a function of C, as in Figure 6, increasing 
S,, shifts the peak deposition higher in the foliage, 
demonstrating a way to determine the shape of deposition 
profile to fit a desired dosage profile. The ratio of velocities 


€ may be chosen to isolate most of the deposited pesticide 
mass within the target subregion, as Figure 7 
demonstrates. In any case, one must calculate the pesticide 
loss by the various possible mechanisms. Figure 8 shows 
the flux of mass loss at the horizontal and vertical 
boundaries, at various values of e, showing a high degree 
of “controllability” of the relative losses through e. Pesticide 
mass loss through mechanisms of vapor production is 
shown in Figure 9, including vapor from the droplet 
spectrum edge, deposited droplets, and free droplets. 








Figure 4—Solution of (7) at three successive times, showing (a) free droplet concentration and (b) deposited droplet concentration. 
Here S, is relatively high. 





Figure 5—Solution of (10) at three successive times, showing (a) free droplet concentration and (b) deposited droplet concentration. 
Here S, is relatively low. 
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Figure 6—For increasing S,, the peak deposition shifts higher into the foliage. 





Figure 7—Deposited droplet concentration at (a) 2,<<U_ and (b) Q.>> U_. 
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Figure 8—Mass loss flux at horizontal (-) and vertical (....) 
boundaries through time: (a) Q2,<<U_, (b) Q,=U_, and (c) 
Q,>>U.. 
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Figure 9—(a) Total vapor mass production from droplet spectrum 
edge r—0 at two successive times, (b) mass deposition at two 
successive times, (c) mass of vapor created from deposited 
droplets at two successive times, and (d) mass of vapor created 
from free droplets at two successive times. 
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Finally, one may investigate the effects of varying the spray 
nozzle, assuming a continuously distributed choice of 
nozzle types based on the initial condition parameter n 
(Figure 10a; see Equation 13). If one has a desired dosage 
profile D(C) as shown by the dashed curve in Figure 10c, 
one may calculate the foliage deposition corresponding to 
each possible n, and compare each resultant deposition 
profile with D using Equation 14. Figure 10c shows two 
nearly optimal profiles, varying n as the control parameter. 


Nozzle Parameter 


0 02 04 06 08 1 


j 


Foliage Deposition 


c 
2 
= 
o 
° 
a 
@ 
ia) 
® 
re) 
S& 
i) 
LL 





Figure 10—(a) Effect of initial conditions variation (due to nozzle 
type) on deposition profile. (b) Resultant dosages for three 
different values of n. (c) Calculated deposition profiles near to 
desired deposition (dashed curve). 
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Discussion 





In this paper we have introduced a preliminary quantitative control parameters and their effects. By disregarding the 
dynamic model for investigating how to optimize the aerial second-order terms in the dynamic equations one may 
application of pesticide to a forest, identifying the various obtain useful, practical solutions. 
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Appendix: Analytical Solutions 





All analytical solutions are calculated by the method of characteristics. 


For the case of only evaporation and deposition, but no movement of the cloud, that is, Q=0, U_=0, e=0 


eS cd (8 
es Fs NAS biog Fist 
49 (fel, (C)e 


which gives the solution 


Only considering streamwise pesticide movement and foliage deposition, that is, Q,=0, e=0, a=0, we have 


OC oc 
pen S reed 8 0 


so 


o(€,6,r) = o(€-U(E)z,C)e 


Thus the deposition rate in the foliage is given by 


oh = k(Q)o(E,S,2) 


or 


g(E,.7) = K(S)f so(E-U(S)0, CO ae’ 


and so the total deposition D(&, C) is 


D(E,¢) = | Saar. 


For the case of no evaporation and no streamwise transport, that is, a=0, U_=0, 


dc Oc 
Beer. k(¢)c =0 
so we have 
e(p,6,t) = (9,6 -Q(p)r)e 
where 








For evaporation with sedimentation, that is, a<0, U_=0, 








dem So eer oc 
— +—| —¢ |+ O(p)— +k(C Jc =0 
$6.49 (o).0(p)E+K(0) 
so 
ah Pp 2 OF, Sig ite ee 
c(p,C, 7) = ——= == ¢) 4) 9° —2at,C - +—2/(p* —2at)2 le 
(p,¢ ) err p a 307 34 (P 
where 
k = 26 +220! 98 3 awl 
rife ca! 0 3a vp ,0 °ed 
8qav2'3 com? e O 
and 
— Qo 
eer 


For the case of saturated conditions, sedimentation, and shearing flow, that is a=0, 


dc dc dc ¥ 
pate ae ele) ap ak\G)e =0 


we obtain on the entire domain 


c(p,€,6,7) = o(p,€ -U_t - A, 6 -Q(p)t)e, 


where 
A= a UE-2(0)2)-U) 
and ; 
1 
p= ap — Q(p)r) ‘ k(¢)] 
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For the general hyperbolic problem, 


oc OCvRre (2 


dc 
ae ) 52+ 35 7 


o}+0(0) S +K(g = 


5 fy PA ee _ 2p" , 20/2 oan)? 
op.8650) = oP of p’ —2at,E+ AF ae Aes 2at) Jereta) 





where 
k = tetl p 3 oe = 
ae ae “e [ee ee ee ale 
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An approach to the problem of modeling atmospheric transport and optimizing 
deposition of pesticide over and through a forest canopy is presented. The 
approach incorporates droplet evaporation, transport of the droplet cloud, 
foliage deposition, and a shearing wind profile. Significant dimensionless 
dynamic parameters are identified, including effective Sherwood numbers, 
Stanton numbers, a distance ratio, and a velocity ratio. Second-order dynamic 
terms in the transport equations are omitted; the resulting hyperbolic 
equations are appropriate for significant limiting conditions of the transport 
parameters, namely infinite Sherwood numbers. Such models prove useful in 
optimization of spray conditions, spray nozzles, and application techniques. 
Analytical solutions are given to the hyperbolic forms of the pesticide transport 
equations. The problem of optimizing the spray application with respect to 
adjustable parameters is described. Various examples are discussed. 
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